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ABSTRACT 


The early twentieth century North American pluvial (1905-1917) was one of the most ex- 
treme wet periods of the last five hundred years and directly led to overly generous water 
allotments in the water-limited American West. Here we examine the causes and dynamics 
of the pluvial event using a combination of observation-based data sets and general circula- 
tion model (GCM) experiments. The character of the moisture surpluses during the pluvial 
differed by region, alternately driven by increased precipitation (the Southwest), low evapo- 
ration from cool temperatures (the Central Plains), or a combination of the two (the Pacific 
Northwest). Cool temperature anomalies covered much of the west and persisted through 
most months, part of a globally extensive period of cooler land and sea surface temperatures 
(SST). Circulation during boreal winter favored increased moisture import and precipita- 
tion in the southwest, while other regions and seasons were characterized by near normal 
or reduced precipitation. Anomalies in the mean circulation, precipitation, and SST fields 
are partially consistent with the relatively weak El Nino forcing during the pluvial, and also 
reflect the impact of positive departures in the Arctic Oscillation that occurred in ten of the 
thirteen pluvial winters. Differences between the reanalysis dataset, an independent statis- 
tical drought model, and GCM simulations highlight some of the remaining uncertainties in 
understanding the full extent of SST forcing of North American hydroclimatic variability. 
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1. Introduction 


The development of western North America (NA) during the twentieth century was largely 
made possible through human appropriation of natural water flows for industrial, municipal, 
and agricultural uses (e.g., Barnett and Pierce 2009; Christensen et al. 2004; Sophocleous 
2010; Reisner 1993; Worster 1992). One set of appropriations is legally formalized under the 
Colorado River Compact (CRC) of 1922 (Christensen et al. 2004; MacDonnell et al. 1995), 
an agreement that apportioned discharge from the Colorado River between the states in the 
Upper (Wyoming, Utah, Colorado, New Mexico) and Lower (California, Arizona, Nevada) 
Colorado River basins (Christensen et al. 2004). The CRC apportionments are based on 
estimated climatological discharge at Lee’s Ferry on the Colorado River of 22 billion cubic 
meters (BCM), using baseline flows from the early twentieth century (Christensen et al. 
2004). As development in the west continued, and as the long term hydroclimate in the west 
was clarified with longer instrumental records and paleoclimate reconstructions (Fye et al. 
2003; Meko et al. 2007; Stockton and Jacoby 1976; Woodhouse et al. 2005), the overly gen- 
erous nature of the original CRC allocations became apparent. For example, mean annual 
discharge at Lee’s Ferry calculated over a much longer interval (1906-2000) was only 18.6 
BCM, ranging in any given year from 6.5 BCM to 29.6 BCM (Christensen et al. 2004). The 
reconstructed discharge extending back to 1512 C.E. is even lower (16.7 BCM), making it 
highly likely that the flows that formed the basis for the CRC were higher than the hydro- 
climatic baseline for the last 500 years (Christensen et al. 2004; Fye et al. 2003; Meko et al. 
2007). 

The exceptionally high flow during the early twentieth century coincided with anomalous 
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wet conditions throughout the west, spanning approximately 1905-1917, a period generally 
referred to as the early twentieth century pluvial (the term pluvial referring to wetter than 
normal conditions) (e.g., Fye et al. 2003, 2004; Woodhouse et al. 2005). This was the most 
persistent pluvial event in the west to occur during the twentieth century, and recent drought 
reconstructions based on networks of tree ring chronologies suggest it may have been the 
wettest period in the west anytime in the last thousand years (Cook et ah 2004). An anal- 
ysis of temperature and precipitation records from the time suggested that the pluvial (as 
reflected in river discharge and drought metrics) arose from a combination of anomalously 
high winter-time precipitation and reduced evaporation from cooler than normal warm sea- 
son temperatures (Woodhouse et ah 2005). 

To date, few studies have discussed the underlying dynamics or causes of the early twen- 
tieth century pluvial. Fye et ah (2004) suggested that anomalously cool temperatures in 
the North Pacific and warm conditions in the tropical Pacific would have favored increased 
moisture flux into the southwest, although this was speculative because of the absence at 
the time of atmospheric circulation datasets covering this time period. Since then, however, 
new datasets and model simulations have become available, leaving us poised for an in-depth 
investigation into the causes of the early twentieth century pluvial in western NA. Here, we 
use available datasets and an ensemble of general circulation model (GCM) simulations to 
investigate the North American pluvial (1905-1917) and assess 1) the relative importance 
of temperature versus precipitation for the pluvial moisture surpluses, 2) the dynamics un- 
derlying these anomalies, and 3) the importance of sea surface temperature (SST) forcing 
during this interval. 
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2. Methods and Data 


To maintain consistency with previous analyses (e.g., Fye et al. 2004), we define the pluvial 
interval as 1905-1917. Over one of our analysis regions (the Central Plains), the pluvial 
moisture surpluses begin several years prior to 1905; the 1905-1917 range, however, is the 
common definition for the pluvial in the literature and this interval includes the major tem- 
poral features of the pluvial at the continental scale. Our analysis will use observation-based 
data sets, proxy based paleo-reconstructions of drought, and a suite of GCM experiments 
to investigate the physics and dynamics underlying the spatial structure and temporal evo- 
lution of the pluvial event. We divide the west into three regions, distinguished (as shown 
below) by the main drivers of their respective pluvial moisture surpluses, their own distinct 
climatologies, and SST-drought teleconnections. Despite these differences, all three regions 
experienced significant wet conditions during the pluvial. The regions are the Southwest 
(SW, 125°W-103°W, 25°N-42°N), the Northwest (NW, 125°W-103°W, 42°N-50°N), and the 
Central Plains (CP, 103°W-90°W, 35°N-50°N). These regions are outlined in Figure 2 and 
other subsequent figures. 

a. Palmer Drought Severity Index 

Droughts and pluvials may be defined in a variety of ways, depending on the research ques- 
tion of interest (Dracup et al. 1980). At the core of all definitions, however, is the concept 
of a moisture deficit (droughts) or surplus (pluvials). While these surpluses and deficits are 
typically viewed primarily as a consequence of moisture supply (i.e. , precipitation), they 
may also strongly depend upon evaporative demand. One drought index that incorporates 


4 


information on moisture supply (via precipitation) and evaporative demand (as a function of 
temperature) is the Palmer Drought Severity Index (PDSI) (Palmer 1965). PDSI is locally 
normalized around a mean of zero with a typical range of -5 to +5; positive values indi- 
cate wetter than normal conditions (pluvials) and negative values indicate drier conditions 
(droughts). PDSI does have some weaknesses. For example, the potential evapotranspira- 
tion calculation embedded within the PDSI uses the Thornthwaite equation (Thornthwaite 
1948), which uses only temperature as an input; in reality, evapotranspiration may also de- 
pend on relative humidity, wind, and cloudiness. Continuous and reliable estimates of these 
other variables can be difficult to obtain, however, and the significance of these other fac- 
tors may be largest in energy limited, rather than moisture limited, regions with relatively 
high background relative humidity (e.g., Hobbins et ah 2008). PDSI also does not, in its 
default form, explicitly model precipitation as snow or accumulated snowpack, which may 
bias results for regions where seasonal snow cover is an important component of the regional 
hydrology. Persistence is built into the PDSI calculation, however, meaning that values over 
one season (e.g., boreal summer) will often partially reflect moisture conditions and climate 
anomalies from previous months and seasons (e.g., St George et al. 2010). Summer (June- 
July-August, JJA) PDSI, the focus of our analysis, should therefore still contain information 
regarding temperature and precipitation anomalies during the antecedent winter and spring, 
in addition to the contemporaneous summer. 

We use two different PDSI datasets in our analysis. The first is the North American 
Drought Atlas (NADA) version 2a (http://www.ncdc.noaa.gov/paleo/pdsi.html) (Cook 
et al. 2007), a tree-ring proxy-based reconstruction of PDSI covering much of North Amer- 
ica. The NADA product reconstructs PDSI for the summer (June-July-August) season 
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using 1,821 tree ring chronologies, over as many as 286 2.5°x2.5° grid-boxes. This product 
is well- validated and versions of the NADA have been used in other studies of NA drought 
variability (Cook et ah 1999, 2004, 2007; Fye et al. 2003; Herweijer et al. 2007). We use the 
NADA to examine the spatial extent and intensity of the pluvial and also place the pluvial 
anomalies within the context of moisture variability over the last five hundred years. 

We also calculate a second set of PD SI values directly from available gridded monthly 
temperature and precipitation data (see below) in order to investigate the relative contribu- 
tion of temperature and precipitation to the PDSI during the pluvial. We set soil moisture 
capacity for the upper and lower layers to 25.4 mm (1 inch) and 127 mm (5 inches) and 
use 1926-2000 as our normalization period (“Climatically Appropriate For Existing Condi- 
tions”). The PDSI calculation is applied to each grid-cell’s temperature and precipitation 
separately; for regional comparisons the new PDSI grids are then spatially averaged over our 
three defined regions. 

b. Temperature and Precipitation 

Gridded temperature and precipitation data are taken from version 2.1 of the Climate Re- 
search Unit (CRU) monthly climate grids (Mitchell and Jones 2005). The CRU data are 
statistically interpolated from monthly station observations to a regular terrestrial grid at 
half degree spatial resolution and monthly temporal resolution, covering the time period 
1901-2002. We use these data to look at seasonal temperature and precipitation anomalies 
during the pluvial and also use them in our own calculation of PDSI. We also use SST data 
from the Hadley Centre (HadlSST; Rayner et al. 2003) and a dataset of global and hemi- 
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spherically averaged temperature for the last 150 years (HadCRUTv3; Brohan et al. 2006). A 
measure of the El Nino Southern Oscillation (ENSO), the NIN03.4 index, is also calculated 
from the HadlSST dataset. Unless otherwise indicated, all temperature and precipitation 
anomalies are expressed relative to a 1961-1990 climatology. We recognize that there may 
be issues using this baseline period because of the strong warming trends over the twentieth 
century. However, we note that 1) this time period is still often used as the standard base- 
line in climate analyses, 2) data during this period is relatively well sampled spatially and 
temporally, 3) persistent drought or pluvial events over NA are generally absent during this 
interval, and 4) it is difficult to develop a comprehensive baseline prior to the pluvial period 
and major warming trends. 

c. Atmospheric Circulation 

Data on atmospheric circulation and dynamics are taken from the Twentieth Century Re- 
analysis Project (20CR) (Cornpo et al. 2006, 2011; Whitaker et al. 2004). This product 
covers the time span 1871-2008, using a data assimilation model forced by observed climate 
forcings and SSTs (from HadlSST). This reanalysis only assimilates surface and sea level 
pressure observations, and versions of the 20CR have already been used to investigate early 
twentieth century circulation features (e.g., Cook et al. 2010; Wood and Overland 2009). 
Because the reanalysis is constructed using an ensemble approach, a useful indicator of un- 
certainty is the ensemble spread in the variables, for example in annual averaged 500 hPa 
heights (Figure 1). The ensemble spread is low where the reanalysis is tightly constrained 
by either a high density of observations (e.g., North America and Europe) or the prescribed 
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SST forcing (e.g., the tropics, where SST forcing of the atmosphere is strong). Over our 
region of interest, the Pacific-North American sector, the spread is on the order of 20-40 
meters over the land, but much higher over the North Pacific. This is not too surprising, 
given the relative paucity of observations available in the North Pacific during this time. 
Despite these caveats, we feel the 20CR is still an appropriate choice for our analysis, given 
the relatively low ensemble spread over NA and our use of the 20CR as supporting analysis 
to our main points. 

d. GCM Experiments 

We also use results from a 16-member ensemble of atmosphere GCM simulations forced 
with observed SSTs to determine the extent to which SST forcing may be able to explain 
climate anomalies during the pluvial. These simulations cover 1856 to the near present, and 
have been previously used to investigate SST forcing of drought over NA with good success 
(Seager et al. 2005b; Seager 2007). The SST forcing comes from Kaplan et al. (1998) for the 
tropical Pacific Ocean for the entire period and, where available, for 1856-1870, and from 
the Hadley center (Rayner et al. 2003) outside of the tropical Pacific from 1871 on. The 
GCM is the Community Climate Model version 3 (CCM3), developed at the National Center 
for Atmospheric Research (Kichl et al. 1998). The model runs are referred to as GOGA for 
Global Ocean Global Atmosphere. All results shown are the mean of a 16-member ensemble, 
representing the SST-forced component of model variability. 


3. Results and Discussion 


a. The Pluvial in a Long-term context 

Anomalous wet conditions were widespread throughout western NA during the pluvial, with 
positive tree-ring reconstructed PD SI anomalies spanning from Mexico to southern Canada 
and from the Pacific coast across the Great Plains (Figure 2). The largest PDSI anomalies 
(+3 and greater) are concentrated along an axis extending from the SW into the NW and 
northern CP regions. Notably, there are no major drought conditions (PDSI <-l) anywhere 
during this time period, at least in the multi-year average. 

All three regions are characterized by high moisture variability and persistent drought 
and pluvial periods over the last 500 years (Figure 3); time series are smoothed (5 year 
lowess filter) to emphasize persistent events. Twentieth-century drought events are well- 
resolved by the PDSI anomalies, including the well-documented droughts in the 1950s (SW) 
and the 1930s ‘Dust Bowl’ (NW and CP), as well as the multi-decadal ‘megadroughts’ in 
previous centuries. Compared to other pluvial intervals over the last five hundred years, the 
early twentieth century event generally stands out through a combination of its intensity, 
duration, and spatial extent. In all three regions, the time evolution of the pluvial indicates 
two wet phases, with a break towards dry or near normal conditions around 1910. The break 
in 1910 corresponds to La Nina conditions in the tropical Pacific, a situation that typically 
suppresses precipitation in SW NA and the southern Great Plains. Over the CP region, 
the pluvial appears to have started earlier than in either the SW or NW, and the post 1910 
phase of the pluvial also appears weaker (drier) in this region. 
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b. Temperature and Precipitation During the Pluvial 


Temperature anomalies were generally cool throughout the west, especially during the spring 
(MAM) and summer (JJA) peak evaporative seasons (Figure 4). The west was also cooler 
than normal during winter (DJF), with the exception of slightly warmer than normal con- 
ditions over California. The largest positive precipitation anomalies occurred in the SW 
during DJF, with increases on the order of 50-60% (Figure 5), and in the SW during MAM, 
although the spring anomalies occur over a much more limited area. There are some minor 
increases over the region of the North American Monsoon, defined here as the areas of Mex- 
ico and the Southwest that receives >50% of their annual precipitation during the boreal 
summer (roughly 112°W-104°W and 22°N-32°N) (Adams and Comrie 1997). In the SW 
regional average, however, these limited positive anomalies during JJA are largely cancelled 
out by the negative anomalies further west. During DJF there were also wet anomalies in 
the CP, but these were relatively low in absolute terms because the annual cycle in precip- 
itation over this region peaks in the summer. What may be even more remarkable is how 
many regions experienced precipitation deficits during the pluvial, especially the SW during 
JJA and autumn (SON), the NW during DJF, and the CP during MAM. This supports the 
hypothesis (Woodhouse et ah 2005) that cool temperature anomalies and low evaporative 
demand may be as important as precipitation for explaining the large moisture surpluses 
reflected in the positive PDSI values. 

A look at the actual temperature (K) and precipitation (mm day" 1 ) anomalies averaged 
over the three regions during the pluvial provide some further insight (Figure 6). Overlain 
in the precipitation plots is a scaled down (60%) version of the NIN03.4 index (dashed 
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line). When NIN03.4 is strongly positive, this is indicative of warm-phase El Nino events 
generally associated with increased winter and spring precipitation in the southwest and 
decreased precipitation during the same seasons in the northwest. During the pluvial, there 
were five significant (>+0.5 standard deviation from the mean) El Nino events: 1905, 1906, 
1912, 1914, and 1915. 

The main evaporative season in all three regions is summer (JJA), and all regions show 
fairly consistent cool anomalies throughout the pluvial during this season (Figure 6, left pan- 
els). In the SW, only 1910 is marginally positive, and the remaining years are all negative, 
with anomalies on the order of -0.5 to -1 K. Absolute anomalies are even cooler in the NW, 
matching or exceeding -1 K in five of the pluvial years. In the CP, the cool anomalies start 
before 1905, averaging about -2 K, and continuing until 1908 with anomalies of about -1 
K, coinciding with the earlier start of the pluvial in this region. Afterwards, temperature 
anomalies are a bit more equivocal, with some major cool years (1912, 1915), but otherwise 
near-normal temperatures. Over the SW and NW, the major precipitation season is DJF; 
over CP annual precipitation peaks during JJA (Figure 6, right panels). Our SW region, 
as defined, does include some areas with significant summer precipitation, especially regions 
associated with the North American Monsoon. Summer precipitation anomalies over the 
SW during the pluvial are muted relative to the anomalies during DJF, however, and we 
also note that we do explicitly account for any summer precipitation changes in our PDSI 
sensitivity analysis (see below). In the SW, there were major positive precipitation anoma- 
lies during both the early and later stages of the pluvial; this contrasts with the NW, which 
showed some minor increases in the beginning, but overall negative precipitation anoma- 
lies throughout. Remarkably, only two of the major precipitation years (1914 and 1915) in 
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the SW actually correspond to El Nino events, despite the significant correlation between 
NIN03.4 and precipitation over the SW during DJF (Pearson’s r=+0.33). Our SW region 
includes some areas where the El Nino/precipitation teleconnections are strongest during 
the spring and autumn (Andrade Jr and Sellers 1988) rather than winter, and we do find 
wet spring (MAM) anomalies in other El Nino years during the pluvial (1905, 1906, 1912), 
although these are muted relative to DJF (not shown). As with the cool temperature anoma- 
lies, the positive precipitation anomalies in the CP begin before 1905. These precipitation 
anomalies persist through the first four years of the pluvial; after 1908 the CP experiences 
precipitation deficits every year except for 1915. 

c. PDSI: Temperature versus Precipitation 

To what extent were the moisture surpluses during the pluvial, as reflected in PDSI, a conse- 
quence of enhanced precipitation versus reduced evaporative demand from cool temperature 
anomalies? To answer this question, we calculate our own PDSI using temperature and 
precipitation data from the CRU climate grids, spatially averaging over the three regions 
(Figure 7, top row). Mean PDSI calculated over the entire pluvial interval is +0.83, +0.73, 
and +0.33 for the SW, NW, and CP, respectively. We expect our calculated PDSI to differ 
somewhat from the values in Figure 3, as the NADA PDSI represent a statistical scaling 
using tree growth as a proxy rather than the direct calculation we conduct here. The cal- 
culated PDSI is generally quite similar to the NADA, showing the two-phase nature of the 
pluvial (pre and post 1910) in all three regions and the early start to the pluvial over the 
CP. 
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To test the importance of temperature versus precipitation as contributing factors to 
summer PDSI during the pluvial, we alternately substitute climatological values (1961-1990) 
instead of observed temperature and precipitation into the PDSI calculation. We substitute 
climatological values for each month (January-December), rather than specific seasons, to 
account for minor but potentially important contributions from anomalies outside the main 
precipitation and evaporation seasons. Substituting climatological temperature (keeping ob- 
served precipitation) into the PDSI calculation results in varying impacts across the three 
regions (Figure 7, center row). The average pluvial PDSI for the SW decreases moderately to 
+0.60, but PDSI during the major pluvial years (1905, 1906, 1907, 1912, 1914, 1915) is only 
slightly diminished. Over the NW, there are major reductions in PDSI across most pluvial 
years, with mean PDSI dropping by over half from +0.73 to +0.31. Even more dramatically, 
mean PDSI switches from positive to negative over the CP (+0.33 to -0.15), suggesting that, 
without the cool temperatures, there would have been drought conditions over this region. 
Calculations with climatological precipitation and observed temperatures (Figure 7, bottom 
row) reduce the original calculated PDSI in the SW by 62% to +0.33. In the NW, mean 
PDSI is reduced to +0.66, a relatively muted response compared to the climatological tem- 
perature scenario, but with the major pluvial years (1907, 1912, 1916) reduced by half. In 
the CP, substitution of climatological precipitation actually enhances the pluvial, increasing 
mean PDSI to +0.84, muting early PDSI anomalies slightly but completely eliminating the 
later drought years during the second half of the pluvial. From these results we conclude 
that the causes of the moisture surpluses varied across these three regions, driven by high 
precipitation (SW), low evaporative demand (CP), or a combination of both (NW). 

Cool temperature anomalies, and the accompanying low evaporative demand, appear 
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to be important for explaining the pluvial moisture surpluses. However, this explanation 
depends on the cool temperatures being an independent factor and not arising from their 
occurrence near the beginning of the twentieth century warming trends or a result of in- 
creased precipitation, which would make things wet and cool by favoring latent over sensible 
heating or increasing cloud cover and shortwave reflection. In the first case, our cool temper- 
atures would simply be a consequence of our relatively warm baseline; if the latter case were 
true, this would mean the temperature anomalies would be a covariate with the precipitation, 
rather than an independent causal factor. A look at global and hemispheric temperatures for 
the late nineteenth to early twentieth century shows that temperature anomalies during the 
pluvial period were cool even relative to previous decades (Figure 8). Thus, an observer in 
1920, looking back, and with no knowledge of the coming warming, would have characterized 
the pluvial decade as anomalously cool compared to previous years. For all three regions 
there is a significant (p<0.05, calculated over 1901-2002) negative relationship between pre- 
cipitation and temperature (Figure 9). However, when the pluvial years are isolated (blue 
dots), we see that temperatures are near normal or cool (left side of the dashed line), regard- 
less of the precipitation anomalies. Because precipitation during the winter and spring may 
be stored as snow pack that carries over into subsequent seasons, the temperature response 
to precipitation may also be lagged by several months. We repeated the correlation analysis, 
this time comparing JJA temperature against antecedent DJF and MAM precipitation. For 
this lagged analysis, only the SW showed a significant correlation (MAM precipitation versus 
JJA temperature: p =-0.23, p<0.05), and with similar results indicating cool temperatures 
during the pluvial that were largely independent of the precipitation anomalies. This gives 
strong evidence to reject the second explanation and conclude that the temperature anoma- 
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lies during the pluvial were largely independent from the precipitation anomalies, allowing 
them to be an independent causal factor for the pluvial moisture surpluses. 

d. Sea Surface Temperatures 

Drought and pluvial events over western NA are largely modulated by variations in SSTs, 
originating primarily from the tropical Pacific (Seager et al. 2005b), part of a zonally and 
hemispherically symmetric pattern of global hydroclimatic variability (Seager et ah 2003, 
2005a). Increased precipitation in the SW is associated with warm-phase El Nino events 
while cold-phase La Nina events typically suppress precipitation over the same region. The 
sign of ENSO-precipitation teleconnections is reversed in the NW, with El Nino events 
leading to drier than normal conditions. The influence of the tropical Pacific is typically 
strongest during boreal winter. 

Composited DJF SST anomalies from all El Nino events (defined as NIN03.4 index > 
+0.5 standard deviation) over the last 130 years are shown in the top panel of Figure 10. 
During a typical El Nino event, warm SST anomalies extend across most of the tropical 
Pacific, flanked by cool SSTs in the extratropical central north and south Pacific ocean 
basins. Warm SSTs also typically occur in a narrow band along the west coast of NA. 
Averaged across all years during the pluvial, SSTs were globally cooler than normal (Figure 
10, bottom panel) and the tropical Pacific is near normal, despite the occurrence of five El 
Nino events (1905, 1906, 1912, 1914, 1915). Even during the pluvial El Nino events, off- 
equatorial SST anomalies in the Pacific deviate from the expected El Nino pattern, especially 
in the extratropical North Pacific which was nearly universally cool across the entire basin 
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(not shown). Coupled with the major precipitation surpluses in the SW during off El Nino 
years, this suggests that El Nino forcing may be insufficient to satisfyingly explain the full 
pluvial anomalies. The North Atlantic was also cooler than normal during the pluvial, a 
condition that is typically associated with increased precipitation in central NA, especially 
during the warm season (Enfield et al. 2001; Kushnir et al. 2010; Mo et ah 2009). 

e. Dynamics 

The impact of various modes of climate variability and atmospheric circulation on climate 
in western NA is well established, including the Pacific North America pattern (PNA) (Wal- 
lace and Gutzler 1981), the Arctic Oscillation (AO) (Hu and Feng 2010; McAfee and Russell 
2008), and the El Nino Southern Oscillation (ENSO) (Andrade Jr and Sellers 1988). What is 
less clear, however, is the extent to which the circulation during the pluvial or other intervals 
represents a forced response to SSTs versus internal atmospheric variability, an important 
distinction as SSTs currently provide the best source of predictability for NA hydroclimate 
(Schubert et ah 2008, 2009). Forcing from tropical Pacific and tropical Atlantic SSTs is 
currently well constrained by numerous studies (e.g., Cook et ah in press; Kushnir et ah 
2010; Seager et ah 2005b) but the role of extratropical SSTs, especially in the extratropical 
Pacific, is still controversial. There is limited empirical and modeling evidence suggesting 
that SSTs in the extratropical Pacific may influence the overlying atmosphere (Frankignoul 
and Sennechael 2007; Kushnir et ah 2002; Peng and Whitaker 1999; Wen et ah 2010). Many 
of the most recent studies suggest that warm SSTs in the Kuroshio extension region may 
force or reinforce an upper level ridge near the Gulf of Alaska during boreal winter, with 
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the atmospheric response typically lagging the ocean by 4-5 months (Peng and Whitaker 
1999; Frankignoul and Sennechael 2007; Wen et al. 2010). Discerning the magnitude and 
significance of the forced response is difficult, however, due to the high internal variability in 
the atmosphere, the fact that extratropical Pacific SST variability is not independent from 
ENSO, and the importance of the atmosphere as the primary driver of SST variability in 
the extratropics (Deser et ah 2010; Junge and Haine 2001; Kushnir et ah 2002). 

Winter (DJF) circulation during the pluvial, as reflected in the 500 hPa geopotential 
height anomalies, is characterized by negative anomalies off the west coast of NA, flanked 
by positive heights over the western extratropical Pacific and over eastern NA (Figure 11, 
top panel). This configuration favors increased moisture advection into the SW along a 
southwesterly track, a pattern highly favorable for so-called “Pineapple Express” winter 
time storm events (Dettinger 2004; Fye et al. 2004). This circulation pattern bears some 
superficial similarities to the canonical circulation response to El Nino forcing, but with some 
substantial differences (Figure 11, middle panel). Specifically, the negative height anomalies 
are weaker and shifted eastward, while the positive heights over NA shift southeast, extend- 
ing over eastern NA and Mexico. The apparent weak circulation response to ENSO during 
the pluvial is consistent with the limited impact of ENSO on the precipitation in the SW 
during DJF and MAM, described previously. The circulation pattern also diverges from the 
expected warm SST-ridge/cold SST-trough response that would be expected given the cold 
North Pacific SSTs (Peng and Whitaker 1999; Frankignoul and Sennechael 2007; Wen et al. 
2010). Low level height anomalies over the subtropical North Atlantic are weak (<5 m) or 
non-existent during MAM and JJA (not shown). The absence of significant high pressure 
anomalies over this region suggests a negligible role for cold North Atlantic SSTs in forcing 
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the precipitation anomalies. 

We calculated the AO index in the 20CR using the leading mode of monthly-mean win- 
tertime sea level pressure following the definition of Thompson and Wallace (1998). The AO 
index from the 20CR correlates well (Spearman Rank Correlation, r=0.98, 1949-2008) with 
the same index calculated from the NCEP-NCAR reanalysis (Kalnay et al. 1996). Spatial 
expression of the AO is similar in both datasets, with positive departures (defined as plus 
one half standard deviation from the mean) in the AO index during boreal winter (DJF) 
characterized by positive height anomalies in the midlatitudes and negative height anoma- 
lies over the Arctic (not shown). Over the Pacific-North America sector (Figure 11, bottom 
panel), this manifests primarily as positive heights over the North Pacific and Eastern NA, 
flanking negative heights over Northwestern NA. Ten years during the thirteen years com- 
prising the pluvial coincided with +AO index departures during DJF (1905, 1906, 1907, 
1908, 1909, 1910, 1911, 1913, 1914, and 1916), and the imprint of these +AO events can 
be clearly seen (Figure 11, top panel). In particular, the positive height anomalies over the 
Western North Pacific and Eastern NA during the pluvial match the location and magnitude 
of the +AO composite anomalies. Temperature and precipitation anomalies associated with 
the +AO composite also show features highly consistent with anomalies during the pluvial 
(Figure 12). Temperatures over Eastern NA are warmer than normal during +AO winters 
and colder over Northeastern Canada and the Northwest; these features are present during 
the pluvial, although the warming is over a more restricted area than would be predicted 
from the AO alone. In addition, +AO winters are associated with modest increases in pre- 
cipitation stretching from the Southwest, across the Central Plains, and into Eastern NA 
just south of the Great Lakes Region. Positive precipitation anomalies occur over this same 
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geographic band during the pluvial. It should be noted that, in a previous study, McAfee 
and Russell (2008) described a minor and largely insignificant relationship between the AO 
and precipitation in the west during boreal winter, apparently at odds with the analysis 
described here. We note, however, that this previous study used a different precipitation 
data set, conducted their analysis over a much shorter temporal interval, and used a differ- 
ent season for their AO index calculation (January- February- March, JFM). Repeating their 
analysis as best we could (JFM AO from 20CR versus DJF precipitation from CRU), we 
were able to largely reproduce their results. Given the methodological differences between 
this study and McAfee and Russell (2008), we feel confident that our own analysis and 
interpretation of the AO impact on precipitation in Western NA is correct. 

/. Modeling 

The AO appears to be an important driver of the pluvial temperature and precipitation 
anomalies during DJF, with El Nino also playing a limited role in the precipitation and 
circulation anomalies during DJF and MAM. To understand the extent to which the plu- 
vial climate anomalies represent a response to external forcing (e.g., SSTs) versus internal 
atmospheric variability, we leverage the GCM simulations described previously to look at 
the climate response over NA to global observed SSTs during the pluvial. Because the AO 
is primarily an internal atmospheric mode of variability, we do not expect that the model 
will be able to reproduce the +AO events and related circulation features during the pluvial. 
However, the model may be able to resolve atmospheric responses to SST forcing originat- 
ing from the extratropical Pacific, as other studies have attempted to do (Frankignoul and 


19 


Sennechael 2007; Kushnir et al. 2002; Peng and Whitaker 1999; Wen et al. 2010). 

Geopotential heights composited from all El Nino years in the GCM ensemble show that 
the model does a good job reproducing the major circulation features associated with El 
Nino in the reanalysis (Figure 13, top panel), although with somewhat weaker amplitudes. 
During the pluvial, however, the height anomaly pattern diverges from the reanalysis and 
seems to simply reflect the imprint of the five El Nino events (Figure 13, bottom panel). 
In the model composite, the positive height anomalies over the western North Pacific and 
eastern NA are absent, and the major region of negative heights is centered over the ocean, 
rather than being shifted to the east over the west coast of NA as it is in the reanalysis. 
Since these features were largely a consequence of the +AO during the pluvial, something 
the model is not capable of reproducing as a response to SST forcing, these differences are 
not surprising. Temperature anomalies from the GCM ensemble (Figure 14) are cool but 
muted compared to the observations and the model also does not reproduce the observed 
spatial structure. The precipitation anomalies from the model largely reflect the influence 
of El Nino events in the mean pluvial circulation, with enhanced precipitation in the SW 
and Mexico in DJF and MAM (Figure 15). These anomalies, however, are muted compared 
to observations (Figure 5), and have slightly different spatial patterns. For example, model 
precipitation anomalies during DJF are shifted too far to the south, and the major pre- 
cipitation increases occur in MAM, rather than DJF. These deviations from the observed 
anomalies are also not surprising in light of the potentially important impact of the +AO 
during the pluvial. As an aside, it should also be noted that two other atmosphere circu- 
lation models, CAM3 and CAM4, are also unable to reproduce the pluvial anomalies (see 
here: http:/ /rainbow. ldeo.colmnbia.edu/%7Ejennic/comparemodel/). 
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In a recent study, Cook et al. (in press) used a statistical drought model to investigate 
SST forcing of persistent multi-year twentieth century hydroclimatic events over western 
NA, including the pluvial. Using statistical modes representing drought variability associ- 
ated with tropical and extratropical Pacific SSTs, they were able to reproduce the magnitude 
and spatial pattern of the pluvial above random noise in over 95% of their model ensemble 
members, with an anomaly correlation between the ensemble median modeled drought pat- 
tern and observations of 0.76. Their conclusion, that knowledge of tropical and extratropical 
Pacific SSTs should be sufficient to predict the pluvial, appears to be at odds with results 
from our SST forced GCM simulations which seem unable to fully resolve the important 
circulation and precipitation anomalies. 

In the statistical model, the mode of drought variability associated with the extratropical 
North Pacific SSTs may simply reflect forcing of SSTs by atmospheric circulation anomalies, 
or could be a high latitude expression of ocean variability linked to dynamics in the tropical 
Pacific (An et al. 2007; Schneider and Cornuellc 2005). This would cause internal variability 
to be mistaken for forcing in the statistical model and overestimate the SST-forced compo- 
nent of the pluvial. On longer timescales than ENSO, Pacific SST variability is dominated 
by the Pacific Decadal Oscillation (PDO), or Pacific Decadal Variability (PDV), which has a 
strong expression in the North Pacific Ocean (Mantua et al. 1997; Zhang et al. 1997). Much 
work has been done, mostly of an observational nature, suggesting that NA hydroclimate 
can be influenced by the PDO (e.g., Gershunov and Barnett 1998; Goodrich 2007; McCabe 
et al. 2004, 2008). For example, differences in a variety of hydroclimatic variables can be seen 
in western NA between positive and negative phases of the PDO, even when only neutral 
ENSO years are considered (e.g., Goodrich 2007). And for various combinations of PDO and 
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ENSO phases, drought and pluvial anomalies can be either amplified or diminished (e.g., 
Kurtzman and Scanlon 2007). However, the PDO is also associated with strong tropical 
Pacific SST anomalies (Zhang et ah 1997) and these PDO-hydroclimate links could be ex- 
plained by the tropical SST anomalies. Indeed, no GCM study to date has demonstrated 
that any appreciable portion of the hydroclimate history of North America is explained as a 
response to extratropical SST anomalies either in the Pacific or Atlantic Oceans. It is worth 
noting, however, that Fye et al. (2004) surmised the atmospheric circulation anomalies dur- 
ing the pluvial using only North Pacific SST information, conjecturing that the combined 
warm tropical Pacific and cold North Pacific would lead to anticyclonic anomalies over the 
western north Pacific and a long-wave trough centered near the west coast of NA. While the 
physical basis for why this would happen is unclear, these circulation features are clearly 
shown in Figure 11. It may be that this GCM, like others, potentially misses an impact 
of SST anomalies in the extratropical Pacific on hydroclimatic variability over NA (which 
are implicitly resolved within a statistical framework) or it could be that the circulation 
anomalies important to the pluvial were a combination of El Nino forcing with significant 
internal atmospheric variability, such as the AO. Only more work, including simulation with 
other GCMs, may be able to resolve this issue. 


4. Conclusions 

Throughout history, persistent periods of extreme climate have significantly impacted the 
functioning of societies, and have often been instrumental in shaping resource use policies and 
societal reorganizations (e.g., Buckley et al. 2010; Hansen and Libecap 2004). One such event, 
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the early twentieth century pluvial, set up unrealistic expectations for water availability 
in western NA, leading to development trajectories that surpassed the long-term support 
capacity defined by the climatology of the region (Christensen et al. 2004). Increasing our 
understanding of the causes and dynamics of this, and other, climate events can help us 
place current and future climate changes in the proper context and inform how we deal with 
these events at the societal level. The specific goal of this study was to investigate the causes 
of the moisture surpluses during the early twentieth century pluvial, and determine how well 
anomalies during that time fit into our understanding of NA hydroclimatic variability. Our 
main results are summarized: 

• Across the west, the origin of the moisture surpluses during the pluvial varied by 
region and can be attributed primarily to increased precipitation (the SW), decreased 
evaporative demand (the CP), or a combination of the two (the NW). 

• El Nino played a partial role in the pluvial circulation anomalies and moisture sur- 
pluses, contributing primarily to increased moisture convergence and precipitation in 
the SW. However, other anomalies in the SST, circulation, and precipitation anomaly 
fields diverged from the expected El Nino response. Forced circulation responses to 
cold extratropical SSTs in the North Pacific were difficult to definitively diagnose either 
within the 20CR or the GOGA model runs. 

• Ten winters during the pluvial were associated with +AO events, and this is reflected 
in much of the circulation, temperature, and precipitation anomaly patterns during 
DJF. The large role played by the AO, an internal mode of atmospheric variability, 
coupled with the relatively weak ENSO and SST forcing at the time, implies that the 
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pluvial was not predictable from SST information. If, however, the AO continues along 
a positive trajectory, as predicted by many models (Miller et al. 2006), it may help 
ameliorate some of the projected drying in the subtropics, including the SW region 
(Seager 2007). 

• However, the intensity and spatial extent of the pluvial can be well reproduced using 
a statistical model with conceptualized tropical and extratropical Pacific SST forcing 
(Cook et al. in press). An independent GCM simulation driven by SST observations 
produces the El Nino response observed during the pluvial, but is incapable of simu- 
lating other important features. 

Studies of pluvial events (e.g., Schubert et al. 2008; Seager et al. 2005b) in the climate 
literature are relatively rare when compared to the wealth of drought investigations. This 
asymmetry is somewhat understandable, given the huge direct and indirect economic costs 
associated with droughts (e.g., Cook et al. 2007). Even climatically advantageous periods 
like the pluvial, however, may have adverse consequences if they are not placed in the proper 
context, as witnessed by the overly generous nature of the CRC (Christensen et al. 2004). 
Extensive research into drought variability over North America has helped illuminate the 
role of SST variability in the ENSO region, and highlighted the importance of La Nina 
events as major drivers of persistent drought in the west (e.g., Seager et al. 2005b). For the 
early twentieth century pluvial, however, our investigation indicates that SSTs in the ENSO 
region had relatively little explanatory power, and highlights the potentially important role of 
internal variability (the AO, cool temperatures, etc) for driving moisture surpluses. Research 
into pluvial dynamics is hampered, however, by the paucity of extended pluvial events that 
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have occurred during the instrumental period. 

The discrepancy between the statistical model and the GCM also highlights some of 
the uncertainties and the often disparate conclusions reached by empirical (McCabe et al. 

2004, 2008) versus model based (Seager et al. 2005b; Seager 2007) investigations of North 
American hydroclimatic variability. Specifically, the two camps disagree on the efficacy of 
extratropical North Pacific forcing of NA hydroclimate with the GCM experiments indicating 
the dominance of tropical forcing. Reducing this key uncertainty will require further studies 
exploiting both empirical analyses and modeling. 
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13 Composite 500 hPa geopotential heights (meters) for all El Nino years (top) 

and all pluvial years (1905-1917; bottom) from the ensemble mean of the 
GCM simulations. 50 

14 Temperature anomalies during 1905-1917 from the ensemble mean of the GCM 

simulations (K), relative to the 1961-1990 ensemble mean. 51 

15 Precipitation anomalies during 1905-1917 from the ensemble mean of the 

GCM simulations, expressed as percent difference relative to the 1961-1990 
ensemble mean. 52 
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20CR Ensemble Spread (m), 500 hPa Hgt (1905-1917) 
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Fig. 1. Ensemble spread (meters) in annual average 500 hPa geopotential heights from the 
20CR during the pluvial interval (1905-1917). 
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NADA v2a PDSI: 1905-1917 
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Fig. 2. Summer season (June- July- August, JJA) tree-ring reconstructed PDSI anomalies 
from the NADA, averaged across all pluvial years (1905-1917). 
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NADA v2a PDSI (JJA): Southwest 



NADA v2a PDSI (JJA): Northwest 
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Fig. 3. Area averaged summer season (June- July- August, JJA) tree-ring reconstructed 
PDSI from the NADA v2a for 1500-2005. All time series are smoothed with a five year 
lowess filter. Time series correspond to the outlined boxes in Figure 1: the Southwest 
(SW; 125°W-103°W, 25°N-42°N), the Northwest (NW; 125°W-103°W, 42°N-50°N), and the 
Central Plains (CP; 103°W-90°W, 35°N-50°N). 
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Fig. 4. Seasonal temperature anomalies during 1905-1917 from the CRU 2.1 climate grids, 
relative to the 1961-1990 mean (K). 
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Fig. 5. Precipitation anomalies during 1905-1917 from the CRU 2.1 climate grids, expressed 
as percent difference relative to the 1961-1990 mean. 
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SW Temp Anomaly (K), JJA 
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Fig. 6. Temperature (K) and precipitation (mm day" 1 ) anomalies during 1905-1917 for the 
three outlined regions. The indicated three month seasons represent the main evaporation 
and precipitation seasons for each region. Dashed lines are a scaled version of the NIN03.4 
index, included to indicate the five El Nino years during the pluvial (1905, 1906, 1912, 1914, 
1915). All temperature and precipitation anomalies are relative to the 1961-1990 mean. 
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SW PDSI (Obs T, Obs P) 
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Fig. 7. Summer season (JJA) PDSI calculated from the CRU 2.1 precipitation and temper- 
ature climate grids averaged over the Southwest (SW, left column), Northwest (NW, center 
column), and Central Plains (CP, right column). For the pluvial period (1905-1917), PDSI is 
alternately calculated with observed temperature and observed precipitation (top row), cli- 
matological temperature and observed precipitation (center row), and observed temperature 
and climatological precipitation (bottom row). 


44 


HadCRU Annual Temp Anomalies (K) 



Fig. 8. Global and Northern Hemisphere averaged surface temperature anomalies (K, from 
HadCRU), 1870-1940. The 1905-1917 pluvial interval (shaded grey) was relatively cool 
compared to either the preceding or following decades. 
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Temp vs Precip: SW, JJA 


Temp vs Precip: NW, JJA 
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Fig. 9. Summer (JJA) temperature versus precipitation anomalies, averaged over the three 
regions for 1901-2002 (from the CRU 2.1 climate grids). Years during the pluvial (1905-1917) 
are colored in blue. Spearman rank correlations between temperature and precipitation are 
shown on each plot. 
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DJF SST Anomaly (K): Pluvial (1905-1917) 
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Fig. 10. Composite DJF sea surface temperature anomalies (K, from HadlSST vl) from all 
El Nino years (n=38) during the instrumental period (top) compared to mean composite sea 
surface temperature anomalies from the thirteen year pluvial interval (1905-1917; bottom). 
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DJF Height Anomaly: 500 hPa, Pluvial (1905-1917) 





Fig. 11. Composite 500 hPa geopotential heights (meters) for all pluvial years (1905-1917; 
top), composite El Nino winters, and composite +AO winters. 
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Fig. 12. Composite temperature and precipitation anomalies for all +AO winters. 
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DJF Height Anomaly (GOGA): 500 hPa, All El Nino Years 



DJF Height Anomaly (GOGA): 500 hPa, Pluvial (1905-1917) 



Fig. 13. Composite 500 hPa geopotential heights (meters) for all El Nino years (top) and 
all pluvial years (1905-1917; bottom) from the ensemble mean of the GCM simulations. 
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Fig. 14. Temperature anomalies during 1905-1917 from the ensemble mean of the GCM 
simulations (K), relative to the 1961-1990 ensemble mean. 
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Fig. 15. Precipitation anomalies during 1905-1917 from the ensemble mean of the GCM 
simulations, expressed as percent difference relative to the 1961-1990 ensemble mean. 
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